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ABSTRACT 



An improved technique for establishing the position, 
dimensions and orientation of a vessel which is essen- 
tially free of internal acoustic discontinuities involves 
scanning the vessel with a beam of ultrasonic energy 
and, from the ultrasonic echoes received in the A-mode 
lines of sight of the ultrasonic beam, establishing a B- 
mode image of the vessel. An angle marker, in the form 
of a short line segment, is established, using the B-mode 
image, on or near a straight axis of the vessel. The 
length of the angle marker is divided into a number of 
sub-sectors, each of which receives approximately the 
same number of A-mode lines of sight of the ultrasonic 
beam. An average A-mode signal is obtained for each 
sub-sector. Comparison of the average A-mode signals 
provides information about the orientation of the angle 
marker relative to the axis of the vessel. Combining the 
average A-mode signals to obtain a composite signal 
and performing edge-detection analysis of the compos- 
ite signal enables the edges of the vessel to be estab- 
lished with precision, and the diameter of the vessel to 
be determined accurately. A single computer can pro- 
grammed to perform the analysis and re-position the 
angle markers (if necessary) on to the axis of the vessel. 
Using information obtained from this technique, accu- 
rate values of the flow rate of a fluid moving within the 
vessel can be obtained when Doppler frequency shift 
measurements are made. 

15 Claims, 4 Drawing Sheets 
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ULTRASONIC METHOD AND APPARAUS FOR 
DETERMINATION OF VESSEL LOCATION AND 
SIZE 

TECHNICAL FIELD 

This invention relates to the measurement of the size 
and orientation of hidden vessels, such as the vessels in 
the human body through which blood flows, using 
ultrasonic echoscopy. It has particular application to 
the measurement of the flow of blood in blood vessels 
using the Doppler frequency shift in conjunction with 
ultrasonic measurement techniques, but it is not limited 
to blood flow measurement. The present invention pro- 
vides a method for improving the current ultrasonic 
techniques for measuring the size and orientation of 
vessels, and hence for obtaining measurements of the 
flow through such vessels. 

BACKGROUND 2 0 

It is now well known that ultrasonic echoscopy tech- 
niques can be used to provide information about an 
object that is not visible to the eye. The basic technique 
of ultrasonic echoscopy involves directing a short pulse 
of ultrasonic energy, typically in the frequency range 25 
from 1 MHz to 30 MHz, into the region of the object 
that is being examined, and observing the energy that is 
reflected, as an echo, from each acoustic impedance 
discontinuity in that region. Each echo received is con 



which is within the vessel being examined. This is 
achieved by fixing the line of sight of the ultrasonic 
transducer and, in the most commonly used version of 
Doppler measurement known as "pulsed Doppler", 
analysing the echoes obtained from the sample volume 
for a limited range of time delays. The Doppler shift in 
the received echoes is averaged in order to calculate the 
average speed of scatterers in the flowing liquid. 

In current applications of the pulsed Doppler tech- 
nique, a small sample volume within the vessel is se- 
lected by the operator of the echoscope, who moves a 
graphical representation of the sample volume over a 
B-mode image of the vessel. In this way, the B-mode 
imaging equipment is used to steer the ultrasonic beam 
and adjust the sample volume delay so that the actual 
sample volume position corresponds to that part of the 
vessel which is to be the subject of the Doppler shift 
measurement. The orientation of the vessel has to be 
known, so that the velocity of the liquid within the 
vessel may be calculated from the well-known Doppler 
equation: 



2 -fa* cos 6 



(1) 



where fo is the Doppler shift frequency, f 0 is the trans- 
mitted frequency, v is the blood velocity, c is the speed 
„ of sound and 9 is the angle between the line of sight of 
verted into an electrical signal and displayed as either a 30 the ultrasonic ^ ^ the direc tion of flow of the 
blip or an intensified spot on a single trace of a cathode 



ray tube or television screen. Such a display of the 
echoes is known as an "A-mode" echograph or echo- 
gram, and is useful in a number of diagnostic techniques 
to locate the boundaries of the object or to provide 35 
other information about the region into which the pulse 
of ultrasonic energy has been directed. 

If a series of adjacent A-mode displays are obtained 
(for example, by physically or electrically moving the 
transmitting transducer which produces the pulses of 40 
ultrasonic energy, or by scanning the direction of trans* 
mission of the pulses of ultrasonic energy), a two-di- 
mensional image of the object under examination may 
be displayed on the cathode ray tube or television 
screen. Such an image or display of acoustic discontinu- 45 
ities is known as a "B-mode" image or display. 

The use of the Doppler frequency shift in the ultra- 
sonic examination of flowing liquids and moving ob- 
jects is also well known. Many echoscopes which per- 
form the B-mode. imaging examination described above 50 
can also perform Doppler frequency shift measure- 
ments in respect of echoes returned from moving ob- 
jects within the region receiving ultrasonic energy from 
the echoscope When the object under examination is a 
blood vessel, measurement of the Doppler shift of ech- 55 
oes from the blood cells within the vessel permits the 
velocity of those blood ceils to be estimated. As pointed 
out by R W Gill, in his article entitled "Measurement of 
Blood Flow by Ultrasound: Accuracy and Sources of 
Error", which was published in Ultrasound in Medicine 60 
and Biology, Volume 11 (1985), pages 625 to 641, it is 
possible to measure the total volume of flow per unit 
time using an ultrasonic examination technique which 
includes the measurement of frequency changes due to 
the Doppler effect. 65 

In ultrasonic examinations including Doppler fre- 
quency shift measurements, it is necessary to obtain 
echoes from a limited volume of the flowing liquid 



liquid. In current implementations of this technique, the 
orientation of the vessel is obtained from observations 
of the graphical representation of the sample volume in 
the B-mode ultrasonic image. 

In the volumetric measurement of flow, a larger sam- 
ple volume is placed to encompass the entire vessel, and 
the total flow is calculated using the relationship: 



flOW 5 



2./ o .cos0 -4 



(2) 



where d is the vessel diameter and is the mean Dop- 
pler shift in frequency. When applying this formula, the 
diameter of the vessel is estimated by the operator, who 
identifies the positions of the two internal vessel walls 
on the B-mode image and places cursors on their im- 
ages. The diameter of the vessel is taken as the distance 
between the cursor positions. This is a difficult measure- 
ment, and because the flow is directly proportional to 
the square of the diameter in the expression for flow, 
errors in the diameter measurement translate into 
greater errors when the flow values are estimated. 

Another factor affecting the accuracy of blood flow 
measurements is the fact that, in humans and animals, 
the diameter of most vessels varies during the cardiac 
cycle. This is particularly so in the case of arteries. 
Hence, for greatest accuracy, the instantaneous values 
of fb and d should be obtained repeatedly and the ex- 
pression in equation (2) should be averaged over several 
cardiac cycles. 

DISCLOSURE OF THE PRESENT INVENTION 

It is an object of the present invention to provide an 
improved ultrasonic echoscopy method and apparatus 
in which the following three pieces of information may 
be obtained in real time: 
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(a) vessel location for accurate placement of the sam- 
ple volume; 

(b) vessel orientation (for the "cos 0" factor correc- 
tion if liquid flow measurements are to be made); 
and 

(c) vessel diameter (for volumetric flow measure- 
ments). 

Preferably, utilising the present invention, the mea- 
surement process can be repeated frequently if any of 
these parameters are likely to change with movement 
(for example, in the cardiac cycle). 

Automatic adjustment of the position of the sample 
volume in ultrasonic echoscopy has been demonstrated 
previously. One technique for doing this is described in 
the paper by J G Davis, K L Richards and D Greene 
entitled "A sample volume tracking unit for pulsed 
Doppler echocardiography", which was published in 
IEEE Transactions on Biomedical Engineering, volume 
BME-26 (1979) pages 285-288. However, that tech- 
nique uses a single line of sight and thus is capable of 
axial shifts only. That and similar prior art techniques 
do not recognise actual structures associated with the 
surrounding tissues, nor do they permit measurement of 
the vessel orientation and diameter. 

The present invention achieves computational effi- 
ciency by using one-dimensional signal processing tech- 
niques. It makes use of two-dimensional data from a 
selected number of adjacent A-mode ultrasonic lines of 
sight in a given B-mode image frame and iterates using 
data from successive images. It uses echoes obtained 
while the echoscope is in its imaging mode, and makes 
use of the facts that (i) images of vessels containing a 
fluid (for example, blood vessels) are normally free of 
internal echoes, and (ii) the vessels on which Doppler 35 
studies are carried out, or in respect of which other 
measurements are likely to be taken, are normally lo- 
cally straight. 

According to the present invention, there is provided 
a method of establishing the dimensions and orientation 
of a vessel which, when subjected to ultrasonic echos- 
copy examination, is found to be essentially free of sub- 
stantial internal acoustic discontinuities, said vessel hav- 
ing an axis which is straight in at least part of the vessel, 
said method comprising the steps of 

(a) obtaining a B-mode ultrasonic echogram image of 
said vessel or the part thereof in which said axis is 
straight from a number of A-mode ultrasonic im- 
ages, each corresponding to a respective A-mode 
line of sight of a beam of ultrasonic energy which is 50 
scanned over said vessel or said part thereof; 

(b) establishing an angle marker in the form of a short 
line segment within or close to the B-mode image 
of said vessel or said part thereof; 

(c) establishing a plurality of adjacent sub-sectors 55 
within the B-mode echogram image of said vessel 
or part thereof; each of said sub-sectors being inter- 
sected by said angle marker; each sub-sector re- 
ceiving substantially the same number of adjacent 
A-mode lines of sight of said scanning beam; 

(d) obtaining an average video-detected A-mode sig- 
nal for each of said sub-sectors, 

(e) performing a comparison of said average A-mode 
signals of the sub-sectors and determining from 
said comparison the difference between the orien- 
tation of the axis of said vessel and said angle 
marker and the orientation of said axis of said ves- 
sel; 



45 



60 



65 



(0 combining said average A-mode signals to obtain a 
composite A-mode signal for said vessel or part 
thereof; and 

(g) performing edge-detection analysis on said com- 
posite A-mode signal, to obtain a value of the diam- 
eter of said vessel or part thereof. 

This method, it will be apparent, may be combined 
with the Doppler measurement technique to obtain an 
accurate measurement of the flow of a liquid within the 
vessel (for example, the flow of blood in an artery, or 
through a cardiac vessel). 

Preferably, the data obtained from the scanning beam 
of ultrasonic energy is obtained in digital form, so that 
the production of the average A-mode signal for each 
sub-sector, the comparison (typically obtained by signal 
superimposition) of such average A-mode signals, the 
difference determinations, the production of a compos- 
ite A-mode signal, and the edge-detection analysis may 
be performed using (a) digital computer, appropriately 
programmed. 

The present invention also encompasses apparatus for 
establishing the orientation and cross-sectional dimen- 
sions of a vessel having an axis, said apparatus compris- 
ing: 

(a) conventional apparatus for generating a beam of 
ultrasonic energy and scanning said beam over at 
least part of said vessel, to obtain a B-mode ultra- 
sonic echogram image of said vessel or part thereof 
from a number of A-mode ultrasonic images, each 
corresponding to a respective A-mode line of sight 
of said beam; 

(b) operator-activated means for establishing an angle 
marker in the form of a short line segment within or 
close to the B-mode image of said vessel or part 
thereof; 

(c) accumulation means for accumulating . echoes 
relating to respective sub-sectors of a plurality of 
sub-sectors of said B-mode image, said sub-sectors 
(i) being adjacent to each other, (ii) being inter- 
sected by said angle marker, and (iii) containing 
substantially equal numbers of adjacent A-mode 
lines of sight of said scanning beam; 

(d) averaging means for producing a respective aver- 
age A-mode signal for each of said sub-sectors; 

(e) comparison means for comparing said average 
A-mode signals and for obtaining differences be- 
tween the average A-mode signals, to enable the 
orientation of said axis relative to said angle marker 
to be computed; 

(0 combining means for combining said average A- 
mode signals to produce a composite A-mode sig- 
nal; and 

(g) analytical means to perform edge analysis on said 
composite A-mode signal and to obtain therefrom a 
value of the diameter of said vessel or part thereof. 

As will be apparent from the comments made above 
concerning the method of the present invention, Dop- 
pler signal-generating and processing equipment may be 
used in conjunction with this apparatus for establishing 
the orientation and cross-sectional dimensions of the 
vessel, to provide apparatus for accurately measuring 
the flow of liquid within the vessel. 

These and other features of the present invention will 
be demonstrated in the following description of an em- 
bodiment of the present invention, which is provided by 
way of example only. In the following description, 
reference will be made to the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS ^.iSS^lZf^^S^ 

FIG. 1 is a block diagram of a flow measurement are parallel to the angle marker line 15), are then accu- 

system incorporating the present invention and a Dop- mulated and an average A-mode video image signal is 

pier frequency shift measurement arrangement, with 5 generated for each sub-sector. In FIG. 2, the average 

the ultrasonic beam scanning arrangement indicated A-mode signal for the sub-sector 23, obtained from 

schematically thereon. echoes received from the region between the lines 27 

FIGS. 2A and 2B illustrate the creation of sub-sectors and 28, is shown as waveform 32. The corresponding 

and contain representations of the average A-mode average A-mode video image signals for the sub-sectors 

signals obtained from the sub-sectors. 10 24, 25 and 26 are shown as signal waveforms 31, 30 and 

FIG. 3 shows waveforms used in part of the edge 29, respectively, 

detection analysis. The waveforms 29, 30, 31 and 32 are displayed in 

FIG. 4 shows another waveform and its processing in FIG. 2B with the axis of the time of arrival of echoes, 

accordance with the preferred edge detection algo- relative to the intersection of the angle marker 15, hori- 

rithm. 15 zontal. It will be seen that although (as expected) the 

FIGS. 5A and 5B are diagrams which illustrate a waveforms 29, 30, 31 and 32 have similar shapes, they 
situation in which the angle marker will be caused to be are displaced horizontally and progressively with re- 
displaced along its length. spect to each other. This apparent offset occurs because 
ncT A TT ™ nccrD TTyrirt vt r\r: toc the marker does not lie parallel to the axis of the 

^t^F^S^SlSSS?^ 20 vessel 13 ' ^ an * le that the ^rresponding points on 

ILLUSTRATED EMBODIMENT the displayed wave forms present relative to the angle 

The apparatus illustrated in FIG. 1 includes a con- marker indicates the angle between the angle marker 

ventional B-mode ultrasonic echograph arrangement and the axis of the vessel. This angle can be used to 

comprising an ultrasonic beam generating, beam trans- determine the orientation of the vessel 13 and/or a new 

mining and echo receiving transducer 10, which is ar- 25 position for the angle marker 15 which is parallel to the 

ranged to scan the ultrasonic beam to and fro through a axis of the vessel 13. 

sector 11. A portion of a tubular vessel 13 (for example, It will also be apparent that each of the waveforms 

an artery) that is to be investigated is located within the 29, 30, 31 and 32 contains low-echogenic regions be- 

sector 11. Echoes received from acoustic discontinuities tween the (parallel) lines 33 and 34, showing that the 

in the path of the lines of sight of the scanned ultrasonic 30 region between the lines 33 and 34 is where the vessel is 

beam are processed in a conventional manner by a time located. The actual position of the intersects between 

gain compensation (TGC) processing unit, non-linear the average A-mode lines of sight signals (assumed to be 

amplification, and video detection stages shown collec- obtained on the mid-line of each sub-sector) and the 

tively as unit 17 in FIG. 1. The B-mode image of the vessel walls is then obtained. Preferably this is done 

region of sector 11 is displayed in a conventional man- 35 using the edge locating algorithm described below with 

ner (on the screen 19 of a television set or cathode ray reference to the composite line of sight for the sector, 

oscilloscope) by the imaging display unit 18. The location of the axis of the vessel 13 relative to the 

A short, linear angle marker is located in the B-mode angle marker is then determined (typically obtained by 

image of the sector 11 in a location corresponding to the fitting a least squares regression straight line, with posi- 

line 15 shown in the sector 11. Thus, effectively, the 40 tion along the angle marker as independent variable and 

operator of the apparatus establishes an angle marker 15 deviation of midline position as dependent variable), 

within a volume 14 of the vessel that is to be investi- This information will enable the angle marker to be 

gated. It will be appreciated that the placing of the repositioned on substantially the axis of the vessel be- 

angle marker assumes knowledge of the position of the fore the next measurement of the vessel 13 is under- 

vessel 13. If that knowledge is not available (for exam- 45 taken. 

pie, if the B-mode image is not entirely clear and the Although the information obtained from the wave- 
vessel 13 cannot be conclusively identified from the forms 29, 30, 31 and 32 is adequate to modify the incli- 
image), the angle marker may be placed by an intelli- nation and range of the angle marker, and may be used 
gent guess of the location of the vessel 13. Using the to estimate the positions of the edges of the vessel 13, 
present invention, it will be possible to progressively 50 the desirability of having a very accurate determination 
re-position the angle marker by an iterative process of vessel size means that usually a more detailed analysis 
until it lies along the axis of the vessel 13. A similar is performed to determine the location of the edges of 
re-positioning of the angle marker will be obtained if the the vessel 13 with precision. 

angle marker is initially alongside or intersecting the A number of algorithms have been established to 

vessel 13. 55 determine accurately the actual positions of the walls of 

Note that the positioning of the angle marker assumes the vessel 13. The preferred algorithm, a one- 

that the vessel 13 is locally straight. dimensional edge detection method, is described below. 

In the illustrated embodiment, as shown in FIG. 2, As noted above, this algorithm is also the one preferred 

the portion of the vessel 13 that is being investigated is when determining the orientation of the vessel, using 

divided into four sub-sectors, referenced 23, 24, 25 and 60 the average A-mode line of sight signals for each sub- 

26. It will be appreciated, however, that any practical sector. 

number of sub-sectors may be adopted. Each sub-sector The average A-mode signals for each of the sub-sec- 

is intersected by the angle marker 15. Each sub-sector is tors (that is, waveforms 29, 30, 31 and 32) are combined 

chosen so that it contains substantially the same number to form a global or composite single A-mode signal, 

of individual A-mode lines of sight of the scanned beam 65 shown as waveform 41 of FIG. 3. 

of ultrasonic energy. With regard to the generation of the combined (or 
The individual A-mode image signals obtained for global, or composite) waveform 41 of FIG. 3, it will be 

each of the sub-sectors, in the region where it is ex- appreciated that if the system illustrated in FIG. 1 is 
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used, it will be necessary to put in a correction factor 
when combining the average A-mode signals, because 
they are not parallel to each other. This is preferably 
done by digitally "re-sampling" the signal waveforms to 
produce the signal waveforms that would be received if 5 
each average A-mode signal was obtained from a beam 
of ultrasonic energy which crosses the vessel at right 
angles to the axis of the vessel. It will also be appreci- 
ated that if, instead of using an ultrasonic transducer 
which scans a beam through a sector, a linear array 10 
transducer is used, so that all of the A-mode lines of 
sight of the scanned beam are parallel to each other, a 
correction will still be required unless the lines of sight 
intersect the vessel axis at right angles. (When produc- 
ing the average A-mode signals for each sub-sector, as 15 
discussed above, using the equipment of FIG. 1, it is not 
necessary to compensate for the change in angle, rela- 
tive to the axis of the vessel, for each line of sight in a 
sub-sector. 

The corrected "raw data" waveform 41 is smoothed 20 
using a Gaussian filter (a suitable degree of smoothing 
corresponds to a Gaussian standard deviation of 1.1 
mm) and the resulting smoothed curve 42 is inspected to 
find its minimum 43. The minimum 43 occurs at the 
location of the axis of the vessel 13. The smoothed 25 
curve 42 is differentiated by the method of first differ- 
ences and the result (waveform or curve 44) is analysed 
to detect a minimum 45, which corresponds to the ves- 
sel wall nearest to the ultrasonic transducer, and a maxi- 
mum 47, which corresponds to the vessel wall farthest 30 
away from the ultrasonic transducer. 

The edge detection algorithm is required to have the 
further capability of distinguishing the vessel wall from 
spurious echoes internal to the vessel, or large specular 
reflectors located outside it. One method of achieving 35 
this is as follows. 

In FIG. 3, using the position corresponding to the 
minimum 43 of the smoothed data curve 42 as a starting 
jgpint, a search is performed in the differentiated data 
curve 44 for minima closer to the transducer than the 40 
starting point, and for maxima further from the trans- 
ducer than the starting point. Within a certain distance 
from the starting point (chosen arbitrarily, but realisti- 
cally, such as 1 2 mm), the largest nearer local minimum 
and largest farther local maximum are found. A second 45 
search^tr!en|locates the respective nearer local minimum^ 
and fanner local maximum nearest the starting point 
and exceeding a suitable fraction (chosen in the light of 
experience and typically about 0.25) of the local mini- 
mum and maximum found earlier. The nearer local 50 
minimum and farther local maximum nearest the start- 
ing point, located by this second search, are chosen as 
the minimum and maximum corresponding the wall 
positions of the vessel 13. 

This selection technique is further illustrated by FIG, 55 
4. In FIG. 4, a plot of the smoothed, differentiated data 
51 is shown, with the minimum of the undifferentiated 
data marked 52. The plot 51 also shows a minimum 
value and three local maxima 53, 54 and 55, which are 
candidates for interpretation as the far vessel wall, are 60 
shown. After measuring the height of the largest local 
maximum 55, a threshold level 56 at a suitable fraction 
(for example, 0.25 of the height of the maximum 55) is 
established. The local maximum 54 nearest the starting . 
point which exceeds this threshold is taken as the wall 65 
position. 

After determining the positions of the walls of the 
vessel 13, the angle marker is normally repositioned to 



8 



lie substantially along the axis of the vessel, and the 
measurements are repeated. 

FIG. 5 illustrates the circumstances in which move- 
ment of the angle marker along its length may be in- 
voked, namely, when the diameter of the vessel, as 
displayed in the image of the vessel, appears to taper 
non-linearly. A linear tapering of the vessel presents no 
problem, but if the image of the vessel tapers non-lin- 
early, this indicates that the image is not entirely along 
the major axis of the vessel. 

In FIG. 5A, the wall positions 61 and 62 of a vessel 
have been determined using an angle marker 63 (shown 
as a dashed line). That the image of the vessel is tapered 
non-linearly is evident from the diameters 64, 65, 66 and 
67, derived from the average A-mode video display 
signals of the sub-sectors 68, 69, 70 and 71, which have 
been plotted as pan of FIG. 5B, If one of the extreme 
diameters 64 or 67 is less than a predetermined fraction 
(such as 0.9) of that predicted by least squares regres- 
sion on the other three diameters, the angle marker is 
moved lengthwise to a new position 72 which is sub- 
stantially on the axis of the vessel, but is farther away 
from the smallest diameter of the vessel, by a suitable 
distance (such as one quarter of its length). The length 
of the angle marker, which defines the length of the 
region of the vessel 13 that is being investigated, will 
normally remain constant. 

From the foregoing description, it will be seen that, in 
its implementation, the present invention records and 
accumulates data from A-mode images obtained using a 
scanning (that is, an imaging) transducer. Preferably the 
data is in digital form. This data is recorded after pass- 
ing through the time gain compensation, non-linear 
amplification and video detection stages of conven- 
tional ultrasound processing. Only A-mode images from 
lines of sight which pass through the angle marker are 
required to be recorded. The part of the scanned beam 
of ultrasonic energy which intersects the angle marker 
is divided into a number of parts (called "sub-sectors", 
typically four in number), each containing an approxi- 
mately equal number of lines of sight of the ultrasonic 
beam. Within each sub-sector, a sum is formed of all of 
the lines of sight within that sub-sector to produce a 
single A-mode video image signal for that sub-sector. In 
accumulating the summation, the time delay of the be- 
ginning of the data from each line of sight is adjusted so 
that the mid-point of each line of sight coincides with 
the angle marker, thus taking into account the assumed 
position and orientation of the vessel. 

The accumulated data segment from each sub-sector 
is analysed using an edge-detection algorithm, to deter- 
mine the location of the near and far vessel walls at the 
centre of the sub-sector. The vessel midline or axis is 
assumed to be halfway between the wall positions deter- 
mined for each sub-sector, measured along the appro- 
priate line of sight (usually the middle line of sight of 
each sub-sector). The measured locations for the several 
midline points (one for each sub-sector) are compared 
with the location of the angle marker and a least squares 
fit is used to calculate the error in range and angulation 
of the vessel relative to the angle marker. 

The data segments for each sub-sector are then added 
together, applying the appropriate correction if the 
lines of sight do not intersect the vessel at right angles to 
its axis, and the wall detection algorithm is applied to 
their sum. The distance between wall positions in this 
last calculation gives the measured vessel diameter. The 
errors in range and angulation (determined using the 



10/24/2004, EAST 



version: 1.4.1 



5,280,787 



least squares regression method mentioned above) are 
converted into range and orientation corrections for the 
angle marker, which are used to establish a new angle 
marker position. The length of the angle marker re- 
mains constant, and its centre normally remains on the 5 
same line of sight. However, if the image of the vessel 
indicates that its diameter tapers non-linearly, the angle 
marker may be moved lengthwise to its new position. 

The new angle marker position is used as input data 
for the next iteration of the vessel tracking algorithm, 10 
which occurs during a subsequent image scan. If the 
initial angle marker position is greatly in error, several 
iterations (one per imaging scan) of the algorithm may 
be necessary to stabilise to the correct value. In a possi- 
ble implementation, the operator places a single cursor 15 
within the vessel of interest, and an angle marker of 
fixed length and orientation (e.g. horizontal) is gener- 
ated as the initial angle marker, which iteratively ad- 
justs itself to the true vessel orientation. 

Those skilled in the art of ultrasonic echoscopy will 20 
recognise that a particular implementation of the inven- 
tion is in conjunction with the known ultrasonic imag- 
ing and pulsed Doppler functions, as shown in FIG. 1. 
These functions may use the same transducer, different 
transducers, or different parts of the same transducer 25 
where that transducer is of an array type. The imaging 
transducer may alternate its lines of sight with the Dop- 
pler lines of sight, or it may assemble a complete image 
during a break in the Doppler data acquisition. 

The vessel orientation information is supplied to the 30 
Doppler module of the machine for substitution in equa- 
tion 1, and so provide correctly scaled velocities of the 
liquid in the vessel 13. In machines which measure volu- 
metric flow, the orientation 6 and diameter d are substi- 
tuted in equation 2 to determine flow. The position of 35 
the sample volume is adjusted for subsequent Doppler 
measurement, and the position and orientation of the 
angle marker become the input values for the next itera- 
tion of the tracking algorithm. 

40 

INDUSTRIAL APPLICABILITY 

Three particularly useful applications of the present 
invention are as follows: 

(a) In studies of the velocity distribution of blood in 
vessels over the cardiac cycle, when the sample 45 
volume is usually smaller than the vessel, this in- 
vention will allow the sample volume to be main- 
tained in a fixed position relative to the (possibly) 
moving vessel, and will automatically calculate cos 

6 for calculating velocity using equation 1. SO 

(b) In volumetric flow studies, when the sample vol- 
ume is usually larger than the vessel, the invention 
will enable the sample volume to be maintained in 
a fixed position relative to the (possibly) moving 
vessel, and will automatically calculate cos 6 and 55 
the vessel diameter d for calculating flow using 
equation 2. 

(c) In Doppler colour flow studies, a two-dimensional 
image is coloured according to the local velocity of 
flow, as described by K. Miyatake, M. Okamoto, 60 
N. Kinoshita, S. Izumi, M. Owa, S. Takao, H. 
Sakakibara and Y. Nimura in their article entitled 
"Clinical applications of a new type of real-time 
two-dimensional flow imaging system", which was 
published in American Journal of Cardiology, vol- 65 
ume 54 (1984), pages 857-868. The present inven- 
tion will allow the use of equation 1 to correct the 
displayed velocity so that the displayed velocity is 
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the actual speed of blood along the vessel, rather 
than the component along the line of sight. 

This list of particularly useful applications of the 
present invention is not intended to be exhaustive. 

Finally, it should be noted that although a specific 
implementation and application of the present invention 
has been illustrated and described above, the present 
invention is not limited to that implementation and ap- 
plication. Variations of and modifications to the present 
invention may be made without departing from the 
present inventive concept. 

We claim: 

1. A method of establishing the dimensions and orien- 
tation of a vessel which, when subjected to ultrasonic 
echoscopy examination, is found to be essentially free of 
substantial internal acoustic discontinuities, said vessel 
having an axis which is straight in least part of the ves- 
sel, said method comprising the steps of 

(a) obtaining a B-mode ultrasonic echogram image, of 
said vessel or the part thereof in which said axis is 
straight from a number of A-mode ultrasonic im- 
ages, each corresponding to a respective A-mode 
line of sight of a beam of ultrasonic energy which is 
scanned over said vessel or said part thereof; 

(b) establishing an angle marker in the form of a short 
line segment within or close to the B-mode image 
of said vessel or said part thereof; 

(c) establishing a plurality of adjacent sub-sectors 
within the B-mode echogram image of said vessel 
or part thereof; each of said sub-sectors being inter- 
sected by said angle marker; each sub-sector re- 
ceiving substantially the same number of adjacent 
A-mode lines of sight of said scanning beam; 

(d) obtaining an average video-detected A-mode sig- 
nal for each of said sub-sectors, 

(e) performing a comparison of said average A-mode 
signals of the sub-sectors and determining from 
said comparison the difference between the orien- 
tation of the axis of said vessel and said angle 
marker and the orientation of said axis of said ves- 
sel; 

(0 combining said average A-mode signals to obtain a 
composite A-mode signal for said vessel or part 
thereof; and 

(g) performing edge-detection analysis on said com- 
posite A-mode signal, to obtain a value of the diam- 
eter of said vessel or part thereof. 
2; A method as defmed in claim 1, in which the step 
of performing edge-detection analysis includes the sub- 
steps of 

(i) smoothing the waveform of the composite A- 
mode signal and identifying a minimum value 
therein which corresponds to the position of the 
axis of said vessel; 

(ii) differentiating the smoothed waveform of the 
composite A-mode signal; 

(iii) identifying a nunimum in the differentiated wave- 
form, which occurs between the generation point 
of the beam of ultrasonic energy and the axis of said 
vessel and which corresponds to one wall of the 
vessel and identifying a maximum in the differenti- 
ated waveform, farther from generation point of 
the beam of ultrasonic energy than the axis of said 
vessel, which corresponds to the other wall of the 
vessel; and 

(iv) adopting the distance between the minimum iden- 
tified in step (iii) and the maximum identified in 
step (iii) as the diameter of the vessel. 
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3. A method as defined in claim 2, including the addi- 
tional step of repeating steps (a) to (g) and using the 
information obtained from step (e) of the last preceding 
sequence of steps (a) to (g) to perform the step of estab- 
lishing the angle marker so that the angle marker is 
positioned substantially on said axis of said vessel. 

4. A method as defined in claim 3, in which, when 
step (e) provides information showing that the B-mode 
image of said vessel indicates that said vessel is non-lin- 10 
early tapered in diameter over at least one of said sub- 
sectors, in the next repeat of said sequence of steps (a) to 
(g), the establishment of the angle marker is effected so 
that the angle marker is effectively moved generally 
lengthwise of its last preceding position. 

5. A method as defined in claim 4, in which said 
B-mode image is displayed on a screen of a television set 
or cathode ray oscilloscope. 

6. A method as defined in claim 4, including the addi- 20 
tional step of performing a Doppler frequency shift 
measurement in respect of fluid flowing through said 
vessel, and using the established orientation and cross- 
sectional dimensions of said vessel to calculate the flow 

25 

rate of said fluid. 

7. A method as defined in claim 1 or claim 2, includ- 
ing the additional step of repeating steps (a) to (g) and 
using the information obtained from step (e) of the last 
preceding sequence of steps (a) to (g) to perform the 30 
step of establishing the angle marker so that the angle 
marker is positioned substantially on said axis of said 
vessel. 

8. A method as defined in claim 7, in which, when 
step (e) provides information showing that the B-mode 
image of said vessel indicates that said vessel is non-lin- 
early tapered in diameter over at least one of said sub- 
sectors, in the next repeat of said sequence of steps (a) to 
(g), the establishment of the angle marker is effected so 40 
that the angle marker is effectively moved generally 
lengthwise of its last preceding position. 

9. A method as defined in claim 1, in which said 
B-mode image is displayed on a screen of a television set 
or cathode ray oscilloscope. 45 

10. A method as defined in claim 1, including the 
additional step of performing a Doppler frequency shift 
measurement in respect of fluid flowing through said 
vessel, and using the established orientation and cross- so 
sectional dimensions of said vessel to calculate the flow 
rate of said fluid. 
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11. Apparatus for establishing the orientation and 
cross-sectional dimensions of a vessel having an axis, 
said apparatus comprising: 

(a) conventional apparatus for generating a beam of 
ultrasonic energy and scanning said beam over at 
least part of said vessel, to obtain a B-mode ultra- 
sonic echogram image of said vessel or part thereof 
from a number of A-mode ultrasonic images, each 
corresponding to a respective A-mode line of sight 
of said beam; 

(b) operator-activated means for establishing an angle 
marker in the form of a short line segment within or 
close to the B-mode image of said vessel or part 
thereof; 

(c) accumulation means for accumulating echoes 
relating to respective sub-sectors of a plurality of 
sub-sectors of said B-mode image, said sub-sectors 
(i) being adjacent to each other, (ii) being inter- 
sected by said angle marker, and (iii) containing 
substantially equal numbers of adjacent A-mode 
lines of sight of said scanning beam; 

(d) averaging means for producing a respective aver- 
age A-mode signal for each of said sub-sectors; 

(e) comparison means for comparing said average 
A-mode signals and for obtaining differences be- 
tween the average A-mode signals, to enable the 
orientation of said axis relative to said angle marker 
to be computed; 

(0 combining means for combining said average A- 
mode signals to produce a composite A-mode sig- 
nal; and 

(g) analytical means to perform edge analysis on said 
composite A-mode signal and to obtain therefrom a 
value of the diameter of said vessel or part thereof. 

12. Apparatus as defined in claim 11, in which said 
accumulation means, averaging means, comparison 
means, combining means and analytical means are con- 
stituted by a single computer programmed to sequen- 
tially perform the functions of said accumulation means, 
averaging means, comparison means, combining means 
and analytical means. 

13. Apparatus as defined in claim 12, including Dop- 
pler frequency shift measuring means for measuring the 
Doppler frequency shift of ultrasonic signals due to the 
flow of liquid through said vessel. 

14. Apparatus as defined in claim 11, including a 
screen for display of said B-mode echogram image. 

15. Apparatus as defined in claim 11, including Dop- 
pler frequency shift measuring means for measuring the 
Doppler frequency shift of ultrasonic signals due to the 
flow of liquid through said vessel. 

***** 
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